Scaling Mobility Patterns and Collective Movements: Deterministic Walks in Lattices 
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Scaling mobility patterns have been widely observed for animals. In this paper, we propose a 
deterministic walk model to understand the scaling mobility patterns, where walkers take the least- 
action walks on a lattice landscape and prey. Scaling laws in the displacement distribution emerge 
when the amount of prey resource approaches the critical point. Around the critical point, our 
model generates ordered collective movements of walkers with a quasi-periodic synchronization of 
walkers' directions. These results indicate that the co-evolution of walkers' least-action behavior and 
the landscape could be a potential origin of not only the individual scaling mobility patterns, but 
also the flocks of animals. Our findings provide a bridge to connect the individual scaling mobility 
patterns and the ordered collective movements. 



PACS numbers: 89.75.Fb, 05.40.Fb, 89.75.Da 
I. INTRODUCTION 

Recently, the scaling properties in mobility patterns 
of animals have attracted increasing attention jl|. The 
traditional scenario about the "nearly random walks" of 
animals is now challenged by the cumulated empirical ob- 
servations, which indicate heavy-tailed displacement dis- 
tributions approximated to a power-law form P{1) ^ Z^". 
Examples include the foraging process and daily move- 
ments of wandering albatrosses 3 j honeybees [3, |j| , spi- 
der monkeys , microzooplanktons @ , marine predators 
0, and so on. These widespread observations imply a 
possible universal mechanism underlying animals' mobil- 
ity patterns. 

Interpretations of the animals' mobility patterns can 
be divided into two classes. One is the optimal search 
strategies [8l-[Tl|. which indicate that animals can max- 
imize the searching efhciency by using power-law move- 
ments. Another is the deterministic walks (DW) (l2| - [l5j 
where a number of preys are randomly distributed on a 
field, and a walker will continuously catch the nearest 
prey from the current position. Recent studies intro- 
duced many real-life factors into the standard framework 
of DW, such as olfactory-driven foraging [l^ and com- 
plex environment effects [iTj . 

To uncover the origin of scaling properties in animal 
mobility, we propose a variant DW model that takes into 
account the regeneration of resources in landscape and 
the least-action movements. Our model can reproduce 
the power-law distribution of displacements, and the scal- 
ing behavior in probability density of having moved a 
certain distance at a certain time, agreeing well with the 
empirical observations. In addition, our model generates 
ordered collective movements of walkers with a quasi- 
periodic synchronization of walkers' directions, indicat- 
ing that the co-evolution of walkers' least-action behavior 
and the landscape could be a potential origin of not only 
the individual mobility patterns, but also the population 



flocks. 



II. MODEL 

The food resources in the real environment can regen- 
erate by themselves as the growth and propagation of 
plants and preys, until reaching a natural limitation of 
abundance. The maximization of foraging benefits and 
minimization of costs (e.g., the least-action movements) 
usually underlie the animals' behavior. Our model takes 
into account these two fundamental ingredients. The en- 
vironment is represented by a two-dimensional Nx N lat- 
tices with non-periodical boundary condition (i.e., walk- 
ers can not go across the boundary), and each lattice 
has prey resource V{i,j) (for the lattice at coordinate 
(«, j)). The maximum prey resource in each lattice is set 
as a fixed value Vm- Different from the standard DW, 
a more realistic case with multi-walkers (the number of 
walkers is denoted by M) as well as their interactions are 
considered in our model. The updating rules about the 
landscape and walkers' positions are as follows: 

(i) At each time step, each walker chooses the nearest 
lattice with the maximum prey resource to occupy at the 
next time step, and if there are more than one possible 
choices, the walker will randomly select one (see Fig. 1). 
The movement is treated as instantaneous, namely the 
diversity of velocity is ignored. The displacement (i.e., 
moving length) I is defined as the geometric distance from 
the current occupied lattice to the next occupied lattice, 
namely, I — ^ (x\ — xqY + (yi ^ 2/o)^i where the coordi- 
nates (a;o, j/o) and {x\^ yi) respectively denote the current 
position and the next position of the walker. 

(ii) The resource V in the current occupied lattice of 
the walker is exhausted by the walker, namely V ^ 0. 

(iii) For each lattice with V < Vm (currently not occu- 
pied by a walker) , the resource increases an unit at each 
time step until V = Vm , representing the regeneration of 
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FIG. 1: (Color online) Illustration of the movement of a 
walker in three successive time steps, where the blue fulled 
circle denotes the current position of the walker, and white 
dashed circles denote all the possible positions in the next 
time step. The values of V in lattices are denoted by different 
gray scales, with pure-white for V = Vm and dark black for 
V ^0. 



prey resources. 

(iv) When the number of walkers M > 1, walkers up- 
date their positions with random order asynchronously 
according to the above algorithms at each time step. 

Because the resource in each lattice regenerates with 
a fixed speed, and each walker consumes at most Vm 
resource at each time step, we define r — MVm/S to ex- 
press the ratio between the total consumption of walkers 
and the total regeneration speed of prey resource in the 
landscape, where S* = x is the area of the land- 
scape. When r = 1, the consumption of prey resource 
is equal to the regeneration speed, and the resource is 
in a critical status. While if r < 1, the system has re- 
dundant resources. In our simulations, M and r are tun- 
able parameters, and the value of Vm is determined by 
Vm = rS/M. 

Noticed that the movements in our model are not 
purely "deterministic" when r < 1 or M > 1. This 
probabilistic property comes from two sources: One is 
the random order in the updating of walkers' positions 
when M > 1, another is that a walker may have more 
than one choices for the next position. Purely determin- 
istic case appears only when r = 1 and M = 1, in which 
the landscape has only one lattice with Vm resource and 
the walker has to periodically repeat its early trajecto- 
ries. In the cases that r — 1 and M > 1, although it 
is possible that walkers exchange their trajectories under 
the treating with random order, the visited time on each 
lattice is still strict periodic if Vm is an integer. 

III. MOVE LENGTH DISTRIBUTION 

In our simulations, the size of landscape is fixed as 
N = 500. We assume the prey resource is full before a 
group of animals come into their habitat, thus we set the 
initial prey resource of each lattice to be Vm , and the ini- 
tial positions of walkers are randomly distributed in the 
landscape. Except for the case of r = 1, our simulations 
show that the move length distribution gets stable after 
the evolution of (i o~r)M steps, so our statistics take 
into account the walkers' movement after ^ o-r)M time 
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FIG. 2: The trajectories of a walker in 5000 consecutive steps 
for different r. Other parameters are AI = 100 and N = 500. 



steps. When r = 1, the number of lattices having maxi- 
mum resource is equal to the number of walkers M in the 
steady state. In this case, each walker has to repeat its 
early trajectory or other walkers' after the first S/M time 
steps (after that time, the consumption equals regener- 
ation), so our statistics take into account the walkers' 
movement after S/M time steps. 

Figure 2 shows the trajectories of a walker for differ- 
ent r, where abundant long-range movements can be ob- 
served when r approaches to 1. The move length distribu- 
tions P{1) for different r and M are respectively shown in 
Fig. 3(a) and 3(b). When r approaches 1, a scaling prop- 
erty of the move length distribution, say P{1) ~ l"", can 
be observed. The analytical result for the move length 
distribution is given in the Appendix, which agrees with 
the simulations. As mentioned above, when r = 1, the 
trajectory is periodic and thus at the last time step of 
a period, the walker returns to its origin, corresponding 
to a generally long displacement with the same order of 
the system size. Therefore, as shown in Fig. 3(a), when 
r = 1, a peak appears at / « iV = 500. The dependence 
between the power-law exponent a and the parameters r 
and M are shown respectively in the two inserts of Fig. 
3(a) and 3(b). Except r = 1, a decreases monotonously 
with the increasing of r. For example, when M = 100, 
a decreases from 3.3 to 2.2 when r changes from 0.80 
to 0.99. This range of a covers almost all the known 
real-world observations of move length distributions of 
animals. This result indicates that the walker is more 
likely to take long-range movement when the prey re- 
source is not rich enough, which is in accordance with 
the experiment on the prey behavior of bumble-bees Q 
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FIG. 3: (Color online) (a) Move length distribution P{1) for 
different r, where the inset shows the dependence between the 
power-law exponent a and r. Other parameters are M = 100 
and N = 500. (b) P{1) for different M, where the inset shows 
the dependence between a and M, and the blue line denotes 
the fitting line with slope —0.192 ± 0.004. Other parameters 
are r = 0.99, and A*' = 500. The red dashed lines in the two 
plots represent a power law with exponent —2. All the data 
points are averaged by 100 independent runs, each of which 
includes 10^ movements. The size of error bars in the two 
insets are smaller than the data points. 



and also is supported by the resent observation on the 
movements patterns of marine predators [18]. Our result 
suggests that the animals living in a habitat with abun- 
dant prey resource will not display scaling property in 
their mobility. 

As shown in the inset of Fig. 3(b) , the relation between 
a and M can be well captured by logarithmic form a ~ 
— InAf. The case oi M — 1 corresponds to a solitary 
animal living in a fixed territory, while M > 1 represents 
the case where several animals share prey resource in the 
same field. This result indicates that the individuals in 
large group are more likely to take long-range movements 
when the resource is not sufficiently rich. 



IV. COLLECTIVE MOVEMENTS 

To our surprise, collective movements are observed 
when M is large and r is close to 1, where walkers may 




FIG. 4: (Color online) A typical ordered marching band of 
walkers generated by our model when t — 50425 {tp — 0.57). 
The red circles and black arrows denote the current positions 
and directions of each walkers respectively. The parameters 
are r = 0.99, M = 500 and N = 500. 
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FIG. 5: The periodic varying of order parameter tp{t) in three 
time periods with parameters r = 0.99, M = 500, and N — 
500. 



line up one or several marching bands, and the walkers in 
the same marching band have similar moving directions 
(approximately perpendicular to the band). A typical 
example is shown in Fig. 4. Sometimes over half walkers 
are in these marching bands. Marching bands are not 
stable. They may suddenly emerge or disappear, may 
grow larger or fall into pieces. 

We define an order parameter iplt) to measure the de- 
gree of synchronization of walkers' directions at time step 
t as the following form (Toj : 



\Em\ 



(1) 



where the velocity vector Vi{t) = {xt+i — Xt)x + {yt+i ~ 
yt)y, and coordinate {x, y) denotes the position of the i- 
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FIG. 6: (Color online) The average value in each period of 
order parameter (ip) (t) for different r, namely (i/j) (t) is the 
average value of tp in the range [t — + ^fj]- Simulations 
run with parameter setting M = 500 and iV — 500. All the 
data points are averaged by 10 independent runs. 
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FIG. 7: (Color online) Comparison of the analytical result 
(Eq. (A5), red line) and the simulation result (black circles) 
when r = 1, M = 1, and A'^ = 500. The term plays a role 
only for large I. 



th walker, x and y respectively denote the unit of velocity 
on the direction x-axis and y-axis. Obviously, ■0 = 1 if 
all walkers have the same direction. 

Figure 5 reports three typical examples of the ip(t) 
curves, respectively for early, middle and relatively later 
stages, where the quasi-periodic behavior with period 
length about S/M can be observed. As shown in Fig. 
6, when r = 0.99, the order parameter will exceed 0.6, 
indicating a strong synchronization of individual's direc- 
tions. Collective movements can be only observed when r 
is close to 1 , which is also the condition for the emergence 
of scaling mobility patterns as mentioned in Section III. 
This property is demonstrated by the simulations shown 
in Fig. 6, where one can see that the steady value of order 
parameter is very sensitive to the parameter r, and when 
r goes down from 1, the collective movements sharply 
disappear. 

Notice that our model does not imply any direct in- 
teraction between walkers. They are driven by the 
co-evolution of resource landscape and the least-action 
movements. This feature is far different from most known 
interpretations on the dynamical mechanisms of animal 
collective movement s [19l-l23| , yet similar to the so-called 
active walk process [24l.l25|. where the macroscopic-level 
structure emerges from the interplay between walkers and 
landscape. The existence of synchronized motions is very 
sensitive to the value of r, suggesting that the food short- 
age may be responsible to the emergence of animal collec- 
tive behaviors, which has also been observed for locusts 
[26j . The existence of the scaling law in displacement dis- 
tribution and the collective movements are, to our sur- 
prise, under the same condition r — > 1. However, these 
two phenomena are not the two sides of a coin, actu- 
ally, they do not straightforwardly depend on each other. 
Whether there exists a certain mechanism underlying the 
co-existence is still an open question to us. 



V. CONCLUSIONS AND DISCUSSIONS 

Our model mimics the mobility patterns of many least- 
action walkers which prey in an landscape with regen- 
erating ability. Scaling laws of move length distribution 
emerges when the regeneration speed of prey resource ap- 
proaches to the critical point that the amount of resource 
is just enough. This result indicates that the mobility 
patterns of animals are sensitive to the environmental 
context (e.g. food resource) jwhich is qualitatively sup- 
ported by real observations [l^ • Our model indicates 
that population also highly affects on the mobility pat- 
terns (see the inset of Fig. 3(b)), which is rarely discussed 
in early DW models. In addition, our model generates 
quasi-periodic collective movements with marching bands 
when r — > 1 indicating that the aggregations of animals 
are more likely to appear with food shortage [1^ [27| . 
which is far different from many known dynamical inter- 
pretations based on the interaction between individuals 
p^l - [23{ . One of the noticeable features in our results is 
the coexistence of both scaling mobility patterns in indi- 
vidual level and collective movement in population level. 
Both phenomena are under almost the same condition: 
the amount of prey resource approaches to the critical 
point, implying that the environment-driven mechanism 
may bridge the scaling individual activity patterns and 
global ordered behaviors. 

Under different parameter settings, our results exhibit 
a wide and gradual spectrum of different mobility fea- 
tures: from the scaling movements to the random-walk- 
like mobility pattern, from ordered collective movements 
to the uncorrelated motions. These results are respec- 
tively well in agreement with different types of real- world 
mobility patterns of animals. However, not all the results 
in our model are fully addressed. Some phenomena, such 
as the logarithmic relation between a and M, and the mi- 
croscopic mechanism in the emergence of such collective 
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movements, arc still open questions to us. 

In a word, although our model is a minimum model 
that many real factors such as the diversity of velocity 
rates and the heterogeneity of environment are not con- 
sidered, the results of our model are generally in agree- 
ment with many wild observations. Our model could be 
helpful in the understanding of the origin of both scaling 
mobility patterns and the ordered collective movements 
of animals. 
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Appendix A: Mean-field Analysis 

When r ^ 1, the analysis can be simplified by the 
periodic movements. At the critical point, the period 
length is equal to S/M and during a period, each lattice 
will be visited exactly once. We introduce a mean-field 
approximation that at any time, the unvisited lattices are 
evenly distributed in the space. After r time steps, the 
number of unvisited lattices is m(T) = S— Mr. Defining 
the normalized move length = l/N = IS~^/'^ where I 
is the real geometric length, the normalized probability 
density of the distance (move length) from an unvisited 
lattice to its nearest unvisited lattice after r steps is: 

p{h,T) « [m(T) - 1] X 2ttU{1 - 7rZ2)™(^)-2. (Al) 



The move length distribution P(Z*) during a period is the 
cumulation of these p{d, r): 

M f^'^ 

Pa*)«^y^ p{U,T)dT. (A2) 

Substitute m(T) = S - Mr into Eq. (Al) and Eq. (A2), 
P(Z*) can be obtained as: 

P{U) « ^[(Ina)-' + (Ina)- V, (A3) 

where o = 1 — ttII. Mostly < 1 and thus ln(l — ttZ^) « 
nil a^nd ~ 1. Therefore P(/*) can be written as 

+ (A4) 

corresponding to the distribution 

p(0~(-r3 + ri). (A5) 



The range of I is limited from 1 to V2S, therefore if S 
is very large, the distribution P{1) is mainly determined 
by the first term Figure 7 reports the analytical 

result Eq. (A5), which agrccis with the simulation well. 
Notice that, the term only show its effect for very 
large I. 
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